Prediction of a novel monoclinic carbon allotrope 
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A novel allotrope of carbon with P2/m symmetry was identified during an ab-initio minima- 
hopping structural search which we call AflO-carbon. This structure is predicted to be more stable 
than graphite at pressures above 14.4 GPa and consists purely of sp^ bonds. It has a high bulk 
modulus and is almost as hard as diamond. A comparison of the simulated X-ray diffraction pattern 
shows a good agreement with experimental results from cold compressed graphite. 



Graphite and diamond are the thermodynamically 
most stable forms of carbon at ambient conditions. How- 
ever, carbon can be found in a vast structural vari- 
ety due to its flexibility to form sp-, sp^- and sp^- 
hybridized bonds: hexagonal diamond, nano-diamond, 
carbon-foams, fullerenes and nanowires are just some ex- 
amples of known carbon allotropes. Experimental evi- 
dences of a novel super-hard carbon phase have been re- 
ported in literature when graphite is compressed at room 
temperatures. Changes in resistivity [H, optical trans- 
mittance 0, Q , optical reflectivity [j| , X-ray diffraction 
(XRD) patterns [11-0] and in the raman spectra Q indi- 
cate a phase transition in the range of 10 to 25 GPa. Re- 
cently, several candidate structures have been proposed 
to match these experimental observations, such as M- 
carbon B, bct-C4-carbon [lo| . VF-carbon (Tlj and Z- 
carbon [8|, |l2| . Although Z-carbon is thermodynamically 
the most promising structure, a final and conclusive de- 
termination has not yet been possible. 

In this article we report on a novel carbon allotrope 
discovered with the recently developed minima hoppin g 
crystal structure prediction method (MHM) [13|, [iJ]. 
The MHM is capable to predict the most stable and 
metastable structures given solely the chemical composi- 
tion of a system. Short molecular dynamics simulations 
are used to escape from local minima, and local geome- 
try relaxations are performed to identify stable configu- 
rations. High efficiency of the escape step is ensured by 
aligning the initial molecular dynamics velocities along 
soft mode direction, and revisiting already explored re- 
gions of the potential energy surface is avoided by a feed- 
back mechanism. The minima hopping method has been 
successfully used in a wide range of applications 



FIG. 1. The structure of MlO-carbon from two different an- 
gles. The left panel shows the 5- and 7-membered rings, while 
the right panel reveals the 6-membered rings. 



projector-augmented wave method was employed as im- 
plemented in the ABINIT code jSc" 
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During our MHM simulations, the energies and 
the Hellman-Feynman forces were evaluated at the 
density functional theory (DFT) level within the lo- 
cal density approximation (LDA), and the all-electron 
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promising structures were further relaxed using norm- 
conserving Hartwigsen-Goedecker-Hutter (HGH) pseu- 
dopotentials . The total energy was converged within 
less than 1 meV per atom by a plane wave cut-off energy 
of 2100 eV and well converged Mankhorst-Pack fc-point 
meshes. We reconfirmed the energy ordering with two 
other exchange-correlation functionals within the gen- 
eralized gradient approximation, namely PBE [23| and 
PBEsol 0. 

Several MHM simulations were carried out with cells 
containing up to 8 atoms at 15 GPa within an uncon- 
strained and thorough structural search, starting from 
different input configurations. We identified a novel, 
monoclinic carbon phase with P2/m symmetry, which 
we call AflO-carbon. It consists solely of sp^ bonds and 
contains 8 atoms per cell. At ambient pressure, the unit 
cell parameters are given by a = 4.080 A, 6 = 2.498 A, 
a = 4.728 A, a = 7 = 90° and /3 = 73.96°. Two 
carbon atoms each occupy the crystallographic 2n sites 
at (-0.1,0.5,-0.113) and (-0.132,0.5,0.421), and the 
2m sites at (-0.333,0,-0.466) and (-0.325,0,-0.117). 
The overall structure is closely related to Af-carbon, also 
consisting of 5- and 7-membered rings along the 6-axis, 
while 6-membered rings are formed along the c-axis. In 
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FIG. 2. Enthalpy differences per atom with respect to 
graphite of several carbon allotropes are shown as a func- 
tion of pressure. The most promising candidate is Z-carbon 
since it crosses the graphite line at roughly 10 GPa. However, 
MlQ-carbon becomes thermodynamically accessible above 
14.4 GPa . 
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FIG. 3. The simulated XRD patterns are compared with the 
experimental results from Ref. 0] of cold compressed graphite. 
While at 13.7 GPa graphite was used to compute the simu- 
lated XRD pattern, MlQ-carbon was used for the simulation 
at 23.9 GPa. 



contrast to M-carbon, the 5-rings share the long edge, 
whereas in MlO-carbon they share the short edge. The 
structure is illustrated in FiglT] 

The dynamical stability was investigated by analyzing 
the phonon dispersion within the whole Brillouin zone. 
The density-functional perturbation theory [25| as im- 
plemented in ABINIT was employed with a 12x12x12 k- 
point sampling and a 4x4x4 g-point mesh. No imaginary 
phonon frequencies were found, confirming the lattice 
stability of the phase. The thermodynamical stability of 
MlO-carbon was investigated by computing its enthalpy 
within a wide pressure range. In Fig. [2] the enthalpies 
of all proposed candidates for cold graphite are plotted 
with respect to graphite as a function of pressure. We ne- 
glected the zero-point vibrational energies in our calcula- 
tions. MlO-carbon becomes enthalpically favorable over 
graphite above a pressure of 14.4 GPa. As expected from 
the structural similarities with M-carbon, both AflO- 
and M-carbon are very close in enthalpy throughout the 
whole pressure range. 

In Table |T] we compare the structural properties of 
AflO-carbon with other carbon allotropes. The bulk 
moduli Bq were computed using the Murnaghan equa- 
tion, and the Vickcr's hardnesses Hy were estimated with 
the method of Gao et al. [2^ . Like all of the investigated 
structures, A/lO-carbon is nearly as hard as diamond and 
has a very high bulk modulus, which could well account 
for ring cracks in diamond anvil cells Furthermore, 
an analysis of the electronic bandstructure was carried 
out, showing that A^flO-carbon is a wide band-gap semi- 
conductor with an indirect DFT gap of 4.4 eV at GPa. 

Finally, we compare the XRD pattern of MlO-carbon 
to experimental measurements from Ref. 0, as illus- 



trated in Fig. [21 The simulated pattern gives a good 
match to the experimental spectrum and could well ex- 
plain the observed changes. It can therefore be ex- 
pected that AflO-carbon is present in samples of cold- 
compressed graphite above 14.4 GPa. However, other 
carbon allotropes show a similarly good agreement and 
thus the XRD pattern alone is by no way a conclusive 
evidence. 

In conclusion we present a novel carbon allotrope 
that we call MlO-carbon. It is a transparent, super- 
hard material which becomes enthalpically favorable over 
graphite at pressures above 14.4 GPa. Both the struc- 
tural and enthalpical properties are very similar to the 



TABLE I. Calculated and experimental data (where available) 
of the bulk moduli Bo (in GPa), Vickers Hardness Hv (in 
GPa) and volumes per atom Vb (in A/atom) at GPa for 
bct-C4, M-, W-, Z-, MlO-carbon and diamond. 



Structure 


Method 


Bo (GPa) 


Hu(GPa.) 


Vo(A) 


bct-C4 


this work 


428.2 


93.5 


5.82 




LDA [llj 


433.7 




5.83 


A/-carbon 


this work 


428.4 


93.9 


5.77 




LDA [9] 


431.2 


83.1 


5.78 


W-carbon 


this work 


427.5 


94.2 


5.75 




LDA [11] 


444.5 




5.76 


Cub-Diamond 


this work 


463.0 


97.8 


5.51 




LDA [llj 


466.3 




5.52 




Expt. 


446^ 


60-120^ 


5.67^ 


Z-carbon 


LDA [8| 


441" 


95.4~ 


5.6" 


A/lO-carbon 


this work 


423.7 


93.5 


5.79 



^ Reference 
^ Reference 
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previously proposed M-carbon [9[. The XRD pattern 
is in good agreement with experimental measurements, 
and although other carbon allotropes are enthalpically 
preferred at lower pressures, it could well be synthesized 
in samples of cold compressed graphite. 
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